The Vertical Structure of Warm Ionised Gas in the Milky Way 



B. M. Gaensler/'^ G. J. Madsen/'^, S. Chatterjee^ and S. A. Mao^ 

^Institute of Astronomy, School of Physics, The University of Sydney, NSW 2006, Australia 
^Department of Astronomy, University of Wisconsin, Madison WI 53706, USA 
'^Harvard-Smithsonian Center for Astrophysics, Cambridge MA 02138, USA 
^ARC Federation Fellow; email: bgaensler@usyd.edu.au 

Abstract: We present a new joint analysis of pulsar dispersion measures and diffuse Ha emission in the Milky 
Way, which we use to derive the density, pressure and filling factor of the thick disk component of the warm ionised 
medium (WIM) as a function of height above the Galactic disk. By excluding sightlines at low Galactic latitude that 
are contaminated by H II regions and spiral arms, we find that the exponential scale-height of free electrons in the 
diffuse WIM is 1830^250 P'-' ^ factor of two larger than has been derived in previous studies. The corresponding 
inconsistent scale heights for dispersion measure and emission measure imply that the vertical profiles of mass 
and pressure in the WIM are decoupled, and that the filling factor of WIM clouds is a geometric response to the 
competing environmental influences of thermal and non-thermal processes. Extrapolating the properties of the 
thick-disk WIM to mid-plane, we infer a volume-averaged electron density 0.014 ± 0.001 cm^"', produced by 
clouds of typical electron density 0.34 ± 0.06 cm~^ with a volume filling factor 0.04 ± 0.01. As one moves off 
the plane, the filling factor increases to a maximum of ~ 30% at a height of ~1-1.5 kpc, before then declining to 
accommodate the increasing presence of hot, coronal gas. Since models for the WIM with a ~ 1 kpc scale-height 
have been widely used to estimate distances to radio pulsars, our revised parameters suggest that the distances to 
many high-latitude pulsars have been substantially underestimated. 
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1 Introduction 

The interstellar medium (ISM) of the Milky Way is a com- 
plex, multi-phase environment. Much of the ISM is an ionised 
gas at temperatures of 10'* — 10*' K, requiring one or more 
vast sources of ongoing energy injection. In the modem pic- 
ture of the ISM (see Ferriere 20011 for a review), ionised gas 
consists of two main phases: a hot "coronal" component at 
temperatures of ~ lO'' — 10*' K, and a warm ionised medium 
(WIM) of temperature ~ 10'' K. 

The WIM is most easily identified in Ha and other op- 
tical recombination lines, which show that much of this gas 
corresponds to H II regions around massive stars. However, 
most H II regions are seen only at low Galactic latitudes. Fur- 
ther from the plane, pulsar dispersion measures (DMs), free- 
free absorption of low-frequency Galactic synchrotron emis- 
sion, interstellar scattering of compact radio sources, and 
faint Ha and H/3 emission all demonstrate the existence of 
a widespread, diffuse WIM, not associated with individual 
stars ( see overviews by iKulkami & Heileslll98l iRevnoldsl 
Il990ah . In the rest of this paper, we use WIM to refer only to 
this diffuse component, disregarding individual H II regions. 

The diffuse WIM has a typical volume-averaged free elec- 
tron density n « 0.01 — .1 cm~'^ and an electron temper- 
ature Te « 8000 K (e.g.. IRevnoldsl [1990^ : IWeisberg et al.l 
I2OO8I) . Scattering and dispersion of high-latitude pulsars demon- 
strate that the diffuse component of the WIM is distributed 
in a thick disk, extending more than a kpc above and be- 
low the G alactic plane with a roughly exponential fall-off 
in density dReadhead & Duffett-SmitlJl975l : lRevnoldsll989l: 
lNordgrenetalJll992h . The power required to maintain the 
ionisation state of the thick-disk WIM is ~ 5x10*^ ergss"^. 



a vast rate of energy injection that can seemingly only be 
supplied by the photo-ionising flux of all the hot stars in 
the Galaxy (Reynolds et al. 1984, 1990b). Questions re- 
main, however, over how these ultraviolet photons propagate 
to large distances ab ove the Galactic plane ( lRevnold3l 1 990d : 
iMiller & Coxlll993l) . and whether the spectrum of the pho- 
ton field produced by massive sta rs is consistent with the 
observed line ratios in the WIM (Reynol ds & TuftellT99^: 
Heiles et al. 1996a; Madsen et al. 2006). Many other edge- 
on spiral galaxies also have a faint, thick WIM di sk, seen in 
deep Ha observations dRandl 1 9981 ; lDettmaJ2(}04h . All these 
results make clear that the WIM is a key part of the feed- 
back process between the stellar and gaseous components of 
a galaxy, and that it connects the energetics and dynamics of 
the disk and the halo. There is thus considerable motivation 
to understand the density, spatial distribution and extent of 
the WIM in the Milky Way and in other galaxies. 

Our own Galaxy is unique in that the DMs of radio pul- 
sars can provide the free-electron column of the WIM for 
thousands of lines of sight, allowing one to develop a smooth 
underlying mod el for the th ree-dimension al Galactic distri - 
bution of n (e.g.. lTavlor"&Co rdes 1993; Go mez et al.l200lh . 
Integrating such a distribution allows one to also predict the 
emission measure (EM) along any sightline. However, even 
when interstellar extinction is taken into account, the pre- 
dicted EMs are vastly smaller than those inferred from dif- 
fuse Ha emission, from low-frequency absorption in the spec- 
tra of non- thermal radio sources, or from free-fre e radio emis- 
sion (e.g.. iPeterson & Webbed I2OO2I ; ISun et aljr200a) . Be- 
cause DMs have a linear dependence on the free electron 
density while EMs depend on the square of the electron den- 
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sity, the inference is that even far from the plane, the WIM is 
not a smooth, homoge neous medium, but is clumped into dis - 
crete clouds or layers ('Revno ldsll977l.ll991bl ; |Pvnzaril993h . 
In the following sections, we review some recent attempts to 
model this structure. 



1.1 Modeling The WIM 
1.1.1 The NE2001 Model 

Recognising the limitations of a smooth, homogeneous WIM, 
Cordes & Lazio (2002, 2003) developed a detailed model, 
termed "NE200r", for the structure of ionised gas in the 
Galaxy. The NE200I model consists of a thin disk of scale 
height 140 pc associated with low-latitude Hll regions; a 
thick-disk WIM layer as discussed above, with a scale height 
of 950 pc; large-scale structure associated with spiral arms; 
cavities and density enhancements corresponding to known 
features in the local ISM; and individual clumps and voids 
needed to model specific regions of enhanced/reduced scat- 
tering or dispersion. The NE2001 model can predict the DM, 
EM, angular broadening, scintillation bandwidth and other 
parameters for any sightline integrated to any distance. It is 
widely used to estimate distances to radio pulsars for which 
there is no other distance indicator other than their DMs. 

The NE200I model attempts to provide a comprehensive 
description of the WIM, and can be continually im proved as 



the Galactic distribution of thermal gas, magnetic fields and 
cosmic rays could then be derived. 

1.1.2 The BMM06 Model 

The arguments in ijl.l.ll suggest that the scale height and 
volume filling factor of the thick-disk WIM need to be re- 
evaluated. A deta i led stu dy of this issue has been carried by 
iBerkhuiisen et al ]( l2006h . hereafter BMM06, who used mea- 
surements of EM and DM for 157 pulsars to derive statistics 
on the density, scale height and filling factor of the WIM. 
BMM06 argued that both the internal cloud electron den- 
sity, A'^ and volume filling factor, /, of the clumpy WIM 
evolve with a source's vertical distance, z, above or below 
the Galactic plane. They fit the data with a model in which 
decays exponentially, while / grows exponentially, with 
z. The expressions for their best fits are: 



N{z) = iVo 
with scale heights 



Hn = TIOI120 pc, Hf = 670Ii3o pc, 
and extrapolated mid-plane values 

No = 0.41 ± 0.06 cm"^ fo = 0.05 ± 0.01. 



(1) 



(2) 



(3) 



aescription ot tne WiJVl, and can be contmually im provea as 

more DMs and other data are obtained. However, as lCordes & Lazid The resulting volume-averaged density of the thick-disk 



l l2002h acknowledge, the NE2001 model fails to correctly 
predict the DMs of some pulsars in high-latitude globular 
clusters. For example, the 22 pulsars in the globular clus- 
ter 47 Tuc have an avera ge DM o f 24.4 pc cm~ ^ and are 
at a d istance of ~ 5 kpc jCamilo et al. 2000; Grat ton et al.l 
l2003h . Pulsars in this cluster are included as an input to 
the NE2001 model, but the prediction of NE2001 is DM = 
42.1 pc cm~^, in si gnificant disagreemen t with the observa- 
tions. More broadlv JLorimer et al.l ( l200d ) have shown that if 
NE200I is applied to the new pulsars subsequently discov- 
ered in the various surveys that used the Parkes multibeam 
receiver, the inferred vertical distribution of these sources in 
the Galaxy has a scale height about a factor of two lower than 
expected. 

The NE2001 model also has other difficulties at high lat- 
itudes. At the north Galactic pole, NE200I predictflEM ^ 
0.05 pc cm"*". Thi s is more than order o f magnitude smaller 
than that observed jRevnoIdsl 1991bl : IPeter son & W ebbeil2002l ; 
iHaffner et alj2003h . suggesting that there is substantial clump- 
ing of the WIM at large distances from the Galactic plane, not 
includ ed in NE200I (see also discussion bv lRevnoldJ 199 IbL 
Il997h . 

Meanwhile, 'Sun et al. ( 2008) have attempted to simul- 
taneously reconcile extragalactic Faraday rotation measure- 
ments and all-sky maps of Galactic radio emission with the 
distribution of n predicted by NE2001. They show that a 
joint fit to all these data requires an anomalously strong halo 
magnetic field, and an unrealistically low scale-height for 
Galactic cosmic rays. Sun et al. ( 2008) argue that if the scale- 
height of the thick-disk WIM in the NE2001 model were in- 
creased from 950 pc to ~ 2 kpc, a self-consistent model for 



WIM, n = fN, remains roughly constant at n ~ 0.02 cm~^ 
up to 2 ~ 1 kpc, implying a much larger effective scale 
height for the WIM layer than had been assumed by all pre- 
vious authors. At higher z, BMM06 argue that the filling 
factor takes on a constant valu e f . 3, wh ile contin- 
ues to decay. As discussed by S un et al.l j2008l) . the NE200I 
model cannot successfully predict the free-free emission seen 
at 23 GHz by the WMAP experiment, but incorporating the 
function f{z) derived by BMM06 (see Eqn. [T] above) pro- 
duces a good match to the data. 

The BMM06 model clearly offers a number of improve- 
ments over previous work in its description of the WIM at 
high z. However, there are several important caveats and 
limitations associated with that study. First, for 95% of the 
pulsars used by BMM06, no independent estimates of dis- 
tance (or hence of z) were available, so BMM06 used the dis- 
tances predicted by the NE200I model, and assumed errors 
of 20% in these distances. However, as we noted in ijl.I.II 
the NE200I model can be in error by much more than this at 
high Galactic latitudes. Conversely, BMM06 explicitly ex- 
cluded from their sample all pulsars in globular clusters or in 
the Magellanic Clouds, eliminating a substantial number of 
sources with reliable, independent distances^ 

A second limitation of BMM06 (and also of the subse- 
quent study by Berkhuijsen & Miiller 2008) is that their data 
sets consist of paired estimates of DM and EM for each pul- 
sar. While the DM is a direct integral of n from the observer 
to the source, the EMs are derived from Hq data, which cor- 
respond to integrals to infinity, attenuated by dust extinction. 



'We have used the software provided at 
|http : / / www .astro ■ Cornell ■ edu/~cordes/NE2001/| 
vl.O. 



subsequent study by IBerkhuiisen & MiilleJ i2008f) has re- 
peated the analysis of BMM06 but only using pulsars with inde- 
pendent distance estimates. H owever, for reasons that they do not 
explain, IBerkhuiisen & Miillerl |2Q0 8) only consider 13 of the 25 
globular clusters Imown to contain radio pulsars, and do not use in- 
formation on pulsars in the Magellanic Clouds. 
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BMM06 scaled down the EM values to account for the fi- 
nite distance of each pulsar, but this calculation propagates 
the significant distance uncertainties in NE2001 into the EM 
values. Furthermore, scaling the EMs requires one to assume 
a WIM scale height, introducing a degeneracy into their sub- 
sequent fitting of the EMs to derive this same scale height. 
BMM06 corrected the EMs for extinction using standard red- 
dening laws, but such corrections assume that dust and gas 
are well mixed along the line of sight. A inhomogeneous 
dust distribution, combined with significant distance uncer- 
tainties for most pulsars, can result in an erroneous extinction 
correction. 

Finally, BMM06 considered pulsars with 60° < £ < 
360° and |6| > 5° (where £ and b are Galactic longitude and 
latitude, respectively). As we will show in ^3.1\ the DMs of 
many pulsars at |6| < 40° are contaminated by Hll regions 
and by other low-latitude structure. When trying to extract 
the properties of the WIM in the thick disk, the DMs and 
EMs of such pulsars will increase the scatter in the data and 
will reduce the signal-to-noise ratio of the underlying signal. 

Overall, we conclude that while the approach of BMM06 
provides several new insights into the behaviour of the WIM 
as a function of z, there are large systematic uncertainties 
that need to be incorporated into the parameters they extract 
from the data. 

1.1.3 Scope Of This Paper 

The discussion in ijl.l.ll and i] 1.1.21 illustrates the strengths 
and weaknesses of the NE2001 and BMM06 models. Specif- 
ically, NE2001 is a detailed description of ionised gas along 
existing sightlines, but has limited predictive power at high 
z where the number of existing sightlines is sparse. BMM06 
provide a comprehensive joint analysis of EMs and DMs, 
and lay out a functional form for the z-dependence of the 
thick-disk WIM that is a significant improvement over that 
of NE2001. However, their analysis suffers from the sig- 
nificant uncertainties introduced by adopting distances from 
NE2001, and by trying to correct EMs for extinction and for 
the finite path length to the pulsars in their sample. 

In this paper, we complement these previous studies, with 
a simple attempt to model the overall structure of the dif- 
fuse WIM. Our main thrust is to re-apply the analysis of 
iRevnoldsl ( Il989l) . who plotted DMj_ = DM. sin|b| vs. z 
for ~ 35 pulsars with known distances, and compared this to 
the predictions of simple geometric models. While this ap- 
proach has been repeated in variou s subseque nt papers (e.g., 
Nordgren et alJl992l:lG6"mez et alJ 2001: Ber khuiisen & Mulled 
2008h . the much larger sample of pulsars now available af- 
fords us the luxury of being able to filter the data to avoid 
contamination by the clumps, spiral arms and H II regions in- 
corporated in NE2001, and to thus boost any signal produced 
by the truly diffuse WIM. In ^ we review the sample of pul- 
sars with known, reliable, distances and define our sample of 
53 sightlines. In Sj3]w6 use a cut on these DMs in Galac- 
tic latitude to derive a substantially revised scale-height of 
the thick-disk WIM, and then combine these DMs with sep- 
arate data on EMs (i.e., not matched to pulsar sightlines, as 
carried out by BMM06 and Berkhuijsen & Miiller 2008) to 
confirm the contention of BMM06 that the data are most sim- 
ply explained if the volume filling factor of this layer evolves 
with z. In Sj4]we discuss the implications of our calculations 
for the distances to pulsars, for the structure of the WIM, 



and for the transition from warm to million-degree gas in the 
Galactic halo. 

2 Pulsars with Reliable Distance Es- 
timates 

Independent distances for pulsars can be determined through 
a variety of techniques. Most commonly, pulsar distances 
are derived kinematically, in which measurements are made 
of the systemic velocities of foreground clouds seen in H I 
absorption (either against emission from the pulsar itself or 
against an associated object such as a supernova remnant). 
A model for Galactic rotation can then be used to convert 
these velocit ies into upper and/or low er limits on the pulsar's 
distance (see lFrail & Weisbergl ll990^for a review). 

While the kinematic approach produces reliable distances 
for many pulsars, there is a key difficulty in using these data 
to study the properties of the thick-disk component of the 
WIM. Since our explicit aim is to study the diffuse compo- 
nent of the WIM, we wish to avoid lines of sight that pass 
through more complicated regions of the ISM. Since kine- 
matic distances rely on the clumpy and complicated nature 
of interstellar gas to produce discrete absorption features, 
sightlines with good Hi absorption measurements are pre- 
cisely those we wish to exclude if we want to study a dif- 
fuse component of the ISM. We thus exclude kinematic dis- 
tances from further consideration, and below focus on three 
remaining samples of pulsars with distances determined in- 
dependently of the foreground ISM properties. In any case, 
even when using only these other distance indicators, we will 
argue in i]3.1.2l that one must exclude lines of sight for which 
|6| < 40° to avoid contamination by Hll regions and spi- 
ral structure. Since H I absorption dis tances are not feasibl e 
for pulsars with |6| > 10° (see, e.g.. lWeisberg et al.ll2008l) . 
this makes moot any concerns as whether H I absorption dis- 
tances should be included in our sample. 

2.1 Pulsars with Parallaxes 

The most reliable pulsar distances come from parallaxes, de- 
rived either through interferometric measur ements or through 
pulsar timing. The pulsar catalogu^^ of iManchester et al.l 
( 2005) lists 34 such pulsars which are also radio emitters (and 
hence have DM measurements). From this sample, we then 
eliminate seven objects for which the uncertainty in the par- 
allax is greater than 33%, and also discard the Vela pulsar 
(PSR B0833^5), because the DM of this nearby pulsar is 
dominated by contributions from its associated young super- 
nova remnant l lBackej|l974lDesai et al.lll992h . 

2.2 Pulsars in Globular Clusters 

Our second reliable sample are radio pulsars in globular clus- 
ters. The on-line compilatiorQ of such sources maintained 
by Paulo Freire lists 137 pulsars in 25 globular clusters. The 

' On-line catalogue at 
|http : / /www. atnf . csiro . au/ research/pulsar /psrcat| 
version 1.33, dated 2008 Jul 19. 

1http:// www . nai c ■ edu/~pf reire/GCpsr . html| , 
version dated 2008 Mar 13 
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DMs of pulsars within each cluster show a small level of vari- 
ation, resulting bo th from ionised gas internal to the cluster 
l lFreire e t al.'200l') and from small-scale angular fluctu ations 



heights z 
30 pc cm^ 



1000 pc, values of DMx saturate at f» 25 — 
correspondingly demonstrating that the electron 



layer falls off in density at these values of z. This saturation 



aye 

in the Galactic DM between adjacent sight lines ( Shishov & Smirnova lue is consistent with the upper limits for the LMC and 



l2002l : lR"ansonil 120071) . However, the fractional variation in 
DM for each cluster is very small: the largest scatter in DM 
is for NGC 6760, for which the fractional variation between 
the two detected pulsars is ~ 2%. We correspondingly con- 
sider each globular cluster as a single datum (regardless of 
how many pulsars the cluster contains), and adopt the mean 
and standard deviation of the pulsar DMs in that cluster as the 
representative value and uncertainty of the DM, respectively, 
for that cluster. Distances to each cluster h ave been take n 
from the on-line versio flof the compilation o f lHarrisl(ll996h . 
The distances were derived from the colour-magnitude dia- 
gram of each cluster by identifying the mean magnitude of 
the horizontal branch population, and as an ensemble should 
be reliable to ~ 15% (C. Heinke, private communication); 
we consequently adopt this as the fractional distance uncer- 
tainty for each cluster. 

2.3 Pulsars in the Magellanic Clouds 

Our final sample ar e pulsars in the Ma gellanic Clouds. The 
pulsar catalogue of [Manchester et al] 1^005) lists 14 radio 
pulsars in the Large Magellanic Cloud (LMC) and five ra- 
dio pulsars in the Sma ll Magellanic Cloud ( SMC) (see also 
[Crawford eta l. 2001 ; M anchester et alj2006l) . Since the SMC 



SMC. In this sense, the data clearly reproduce the find ings of 
many previous studies (Revnolds 1989, 1991a; Savage et al.l 
1990; Bhattacharya & Verbunt 1991; Nordgren et al. 199? 
Mill er & Cox. 19 93; Gomez et al. 2001, ;,Cordes & Lazia20^ ; 
Hill et alj2005l; lBerkhuiisen & Mulle fcoOSL BMM06). 

The overall distribution of ionised gas in the thick disk 
is reasonably described by a single plane-parallel layer (e.g., 
Haffner et al. 2003). To confirm that our data are consistent 
with such a geometry, and to estimate the scale height of this 
structure, we compare our measurements to a planar distribu- 
tion of free electron density, n{z), which falls off exponen- 
tially with increasing z: 



n{z) = no e 



(5) 



where no is the extrapolated density of ionised gas at mid- 
plane and H„ is the corresponding scale height. Integrating 
Equation (O out to a given distance then gives the predicted 
distribution of DM_l (z): 



DM±{z) = no H„ 



(l-e-*) 



(6) 



3.1.1 Exponential Fit To AU Data 



and LMC themselves have significant DM contributions dGaensler e^ al.^e fit the 5 1 measurements in Figur efflto Equation 

|2005; Manchester etal. 2006; Mao etal. 2008), these sys- ^ — ■■ — ^ 

tems can only provide an upper limit on the Galactic DM 



contribution along these sight lines. For the LMC, the low- 
est DM is that of PSR J0451-67, but there is uncertainty as 
to whether this pulsar i s actually in the LMC, o r is a Galac- 
tic foreground source ( [Manchester et al.ll2006l) . We corre- 
spondingly adopt the next highest DM, 65.8 pc cm~^ for 
PSR J0449-703I, as our upper limit in the direction of the 
LMC. For the SMC, the lowest DM is 70 pc cm"^ for PSR 
J0045-7042. We assume dis tances to the LMC and SM C of 
50 and 61 kpc, respectively jWalkei[l2003l ; [Alve j[2004l . and 
references therein). 



3 Analysis 

3.1 The Distribution of Dispersion Mea- 
sures 



The DM of a pulsar is defined as 



DM 



n{l) dl, 



(4) 



where d is the distance to the pulsar and n(l) is the free elec- 
tron density along a line element dl. The distributions of 
DMx = DM. sin |6| as a function of z for the 51 measure- 
ments ( i]2.1l & i]2.2b and two upper limits ( i]2.3l l are shown 
in Figure [T] It is immediately clear that for pulsars with 
vertical heights in the range IQ i z 1000 pc, the av- 
erage value of DMx steadily increases with z, indicating 
that such sources sit within the Galactic electron layer. For 



jhtt p : / /phy swww . mcmaster . ca/^harris/mwgc . dat] , 
version dated 2003 Feb 



using a Levenberg-Marquardt algorithm (iPress et al.llI992 
(the two upper limits from the Magellanic Clouds are ex- 
cluded from the fit, but are used as independent tests on the 
validity of any fitted model). Note that since the errors in 
DMx are negligible compared to the errors in z, we neglect 
the former and adopt DMx as our independent variable, with 
each datum weighted by the inverse square of the error in z. 

A fit of all the data to Equation l|6) does not converge. 
This is due to the very small errors in z for pulsars with dis- 
tances from parallaxes: since there is no single curve that 
passes through all these points, the resulting is always 
high, with no set of fit parameters that correspond to a signif- 
icant minimum. This problem can be mitigated by recognis- 
ing that although the errors in distances to individual pulsars 
are indeed small, there is an additional systematic effect, in 
that random clumps of electrons along individual sightlines 
cause each measurement to deviate from t he smooth u nder- 
lying model ( Gomez et al. 2001 ; Cordes & Lazid[2002[) . 

Rather than explicitly incorporate individual clumps to- 
ward each pulsar as for the NE2001 model, we represent 
the overall systematic limitations of the data by including in 
the fit an additional fractional uncertainty (added in quadra- 
ture) to the value of z for each pulsar (see, e.g.. Savage et al] 
'1990). Adding successively larger fractional errors estab- 
lishes that an additional 10% systematic uncertainty allows 
the fit to Equation ^ to converge, with a reduced chi-squared, 
Xr, of 12.3 for 49 degrees of freedom. The best-fit parame- 
ters are //„ = pc and no = O.OSlti^ gga cm-^ 
where the uncertainties have been determined via a bootstrap 
algorithm (Efron & Tibshirani 1993), and correspond to the 
most compact 68% interval of bootstrap realisations around 
the best-fit values. The corresponding curve is plotted as a 
dashed line in Figure [H this fit is consistent with the up- 
per limits on DMx for the Magellanic Clouds, and is quite 
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similar to the tiiick-disk component of NE2001, for whicli 
H„ = 950 pc and no — 0.035 cm~^ (plotted as a dot-dashed 
line in Fig.CjQ 

3.1.2 Exponential Fit To ffigh-Latitude Data Only 

While the results in i]3. 1. ll demonstrate that we can reproduce 
the findings of previous studies, the quality of the overall fit 
is reasonably poor. In particular, the best-fit curve substan- 
tially over-estimates DMx for most pulsars at 2; > 1000 pc. 
A possible cause of this discrepancy is indicated in Figure[2] 
where we plot the distribution of our sample on the sky in 
Galactic coordinates. This demonstrates that there is a sig- 
nificant concentration of sources (mainly pulsars in globular 
clusters) in the Galactic plane and toward the Galactic Cen- 
tre, for which the DMs are relatively high. This is simply 
understood as resulting from these pulsars being viewed not 
just through the thick-disk component of the WIM, but also 
through an additional thin-disk layer, corresponding to the 
dense, turbulent, ionised gas in individual Hll regions and 
supershells. The greyscale in Figure [2] clearly demonstrates 
this effect (see also Fig. 6 of BMM06), showing that low- 
latitude pulsars usually lie behind complicated regions of en- 
hanced Ha and EM. 

We have correspondingly colour-coded the data in Fig- 
ure[T]by Galactic latitude. It can indeed be seen that the red 
and orange points tend to lie above the blue and green ones, 
demonstrating that for a given z, sources at low |6| tend to 
have higher values of DMx than at high |6j. To isolate the 
thick disk component, we thus fit to successive subsets of 
the data, selected by adopting a minimum threshold for 
|6|min, and thereby excluding sources contaminated by the 
thin-disk component of the WIM. 

The results of these fits for 0° < |&|min < 75° are shown 
in Figure[3] This plot demonstrates that the quality of the fit 
is relatively poor, Xr ~ 12 — 13, for 0° < |6|miii < 30°, 
but is significantly improved, Xr ~ 4 — 5, for fits with 
40° < |6|miii < 65° (above this value, the sample size be- 
comes too small to be meaningfully analysed). We adopt 
as our best description of ionised gas in the thick disk the 
lower latitude threshold that maximises the size of the sam- 
ple while maintaining a good fit. This value is |6|min = 40°, 
for which Xr = 5.2 (for 13 degrees of freedom); the corre- 
sponding 15 systems are indicated as blue data points in Fig- 
uresEandU The fit to these data yields H„ = 1830l25o 
and no = 0.014 ± 0.001 cm~^, shown by the solid line 
in Figure [1] The bootstrap approach incorporates the added 
uncertainties due to systematic deviations from the smooth 
underlying model. Thus although Xr 3> 1, the errors in Hn 
and no meaningfully characterise the confidence intervals on 
these parameters. 

Our estimate of Hn is substantially larger than that pro- 
posed by previous authors, as summarised in Table [T] How- 
ever, our is clearly a good description of virtually all the 



ICordes & Lazid j2002h adopt a distribution n{z) oc 
sech^ (z/Hn) rather than the exponential form of Eqn. (Sj- 
but the two functions are viitually indistinguishable, especially at 
high z. 

'We note that the best-fit parameters do not change appreciably 
for larger values of |fe|min- f™ example, for Iblmin = 65°, the 
fitted parameters are H„ = ISTOljso 1"^ "0 = 0.0139 ± 
0.001 cm-3, with Xr = 6.0. 



high latitude data, especially those pulsars at high z. The im- 
plied integrated electron column density toward the Galactic 
poles is DMo = noH„ = 25.97^2 35 P'- cm~^, consistent 
with the upper limits from pulsars in the Magellanic Clouds, 
and comparable to the integrated vertical DM found in most 
previous studies. 

The reason for our much larger estimate of H„ is clear 
from Figure[T] the high-z pulsars and Magellanic Cloud data 
together constrain DMo for any model. Curves with high 
values of no must thus turn over early to match this required 
value for DMq, while those with low no can extend to high 
values of z before saturating. When fitting all data (shown 
by the dashed line), the group of red and orange points at 
300 pc ^ z i 1000 pc forces the best fit to have a com- 
paratively high value of no, and hence a low value for Hn. 
The blue data points all systematically lie below the other 
measurements, and the corresponding fit (shown by the solid 
line) must then have a large value of Hn to converge to the 
required value for DMq . 

3.2 Scale Height of Emission Measure 

In i]3.1l we focused exclusively on pulsar DMs. However, an 
independent measure of the free electron column along the 
line of sight comes from measurements of the Ha surface 
brightness, which can then be used to infer the corresponding 
EM, 

poo 

EM= / n^{l)dl. (7) 
Jo 

The integral in Equation l|7j differs from that for pulsar DMs 
in Equation ^ in that the former does not terminate at any 
specific distance, and that the Ha emission must be corrected 
for dust extinction to obtain the corresponding EM. Thus 
EMs and DMs along the same sightlines can only be mean- 
ingfully compared for sources for which both \b\ and z are 
simultaneously high, so that extinction is low and the sources 
are entirely above the WIM layerO This means that EMs to- 
ward pulsars ca n be used to infer volume filling fractions of 
ionised gas (see lRevno idi ll99Ibh . but cannot be used to in- 
dependently estimate the gas scale height, as was done for 
DMs in Figure [T] 

A separate approach to the EM scale height has been en- 
abled by the wide-field imaging spectro scopy of the WIM by 
the Wisconsin Ha Mapper (WHAM). iHaffner et all ( Il99a) 
used WHAM to map the intensity of diffuse Ha (i.e., that 
emission not associated with individual H II regions) for the 
restricted velocity range corresponding to the Perseus spiral 
arm as a function of b, and found that this intensity dropped 
exponentially by a factor of 100 between b = —15° and 
b — —35°. For a planar distribution of ionised gas with an 
exponential density profile: 

n'(2) = (n')oe"^, (8) 

and assuming a distance to the Perseus arm of 2.5 kpc. lHaffner et al.l 
( 1999) derived a scale height for the EM layer of H„2 = 
500 ± 40 pc. The distance to the Perseus arm has sub- 
sequentl y been determi ned by trigonometric parallax to be 
2.0 kpc Ixu et al.ll2006h . so we revise the published value of 
the EM scale height in Perseus to _H'„2 = 400 ± 30 pc. 

'^As discussed in ^1.1.21 BMM06 and iBerkhuiisen & MuUeil 
12008) attempted to correct for these effects for sources at lower |b| 
and z. 
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3.3 The Filling Factor of Ionised Gas 

We have above derived two independent estimates of the WIM 
scale height: Hn — I830I250 P*- using the dispersion mea- 
sures of pulsars at known distances in i ]3.1.2l compared to 
H^2 — 400 ± 30 pc from the latitude-distribution of Ha 
emission in and above the Perseus arm, as discussed in i)3.2l 
Comparison of Equations ^ and {Sjl clearly requires that 
H„/H^2 = 2, which is in strong disagreement with the ob- 
served value H„/H„2 — 4.7lQ g. As discussed in i[T] such 
discrepancies can be resolved by introducing the volume fill- 
ing factor, /, such that along a given sightline, a fraction / of 
the line of sight intersects clouds of uniform electron density 
N, with the rest of the path not contributing to the observed 
EM or DM (e.g., because it is occupied by neutral gas)Q In 
some local region, we then have: 



n{z) = f{z) N{z), 



(9) 



where TV is the electron density within a WIM cloud, and n is 
the mean electron density, averaged over some scale signifi- 
cantly larger than the av erage cloud sep aration. As originally 
proposed by Kulkarni & HeilesI | |19871 |1988) and discussed 
extensively by BMM06, we can now separately consider in- 
dependent distributions of / and A'^, both as a function of z: 



f{z)=foe'^"f, 
N{z) =Noe"^ 



(10) 
(11) 



where /o and TVo are the extrapolated filling factor and cloud 
density at mid-plane, respectively, and Hf and Hn are the 
corresponding scale heights of the two distributions. Note 
that / is an increasing exponential function of 2 in Equa- 
tion diet , since we expect the WIM to increasingly dom- 
inate over neutral gas at progressiv ely higher values of z 
jKulkami & Heileslll987l: |de Avillez 2000). 

Substituting Equations JlOl l and i ll 11 1 into Equation ([9]l, 
we then obtain: 



n{z) = /o 
so that from Equation {5) we then infer: 



1 



no = 
1 

Hn 



foNo, 
1 



Similarly, for density squared, we find: 

n\z)^foNSe~i['^~ 
which implies: 



{n )o = fo No, 
1 _ 2 1 
Hn2 Hn Hf 



(12) 

(13) 
(14) 

(15) 

(16) 

(17) 



'We have defined a line-of-sight filling factor, which is not nec- 
essarily equivalent to the three-dimensional volume filling factor. 
However, for reasonable assumptions about a clumpy distributed 
medium, the two different filling factors are equivalent (see §2.1 of 
BMM06). 



Equations l ll4b and J17b can then be solved jointly to 
yield: 



J/„2 H„ H„2 Hn 

tlN ~ — TZ ,-H/ — 



Hn — Hn' 



Hn — 2J/„2 



(18) 



Our adopted values Hn = 1830^250 and H„2 = 400 ± 
30 pc then imply separate scale heights. 



Hn = 5WlLtl pc, 
Hf = 700t7SJ'' pc. 



(19) 
(20) 



Note that despite the broadly comparable values of Hn and 
Hf, the error bars on these quantities are not gaussian, and 
we can exclude the possibility that Hn ~ Hf at high confi- 
dence. 

We can separately determine the normalisation factors fo 
and A^o, as defined in Equations ilO\ and dl It . respectively, 
as follows. We determined in i|3.1.2l that pulsar data im- 
ply an extrapolatecPI mid-plane mean electron density no — 
0.014 ± 0.001 cm"''. The mean vertical component of the 
Ha intensity at high latitude, assuming a planar geometry fo r 

the emitting region, is 0.68 rayleighs jHaffner et al.ll200"3h. 

For electron temperatures in the range 8000 to 12 000 K dRevnolds et al.l 
[1999; Haff ner et alj|l99^ , this corresponds to a total verti- 
cal emission measure EMo ~ 1.9 ± 0.3 pc cm^'' (the cor- 
rection for dust extinction at high latitudes is a t the level of 
~ 2% and can be ignored; ISchlegel et alj|l99^ . If we com- 
bine Equations l|7j, ^ & l ll6l l. we determine: 



EMo = /oA^c?i^„2. 



(21) 



Assuming that the observed Ha emis sion at high lati tudes 
originates from the diffuse WIM (see iHill et"aiTl2008ll . we 
can combine Equations l IlBt & | |21| | to yield: 



A^o 



fo 



fom ^ EMo 
foNo Hn2 no 
ifoNof 



0.34 ± 0.06 cm" 



no 



H„ 



EMo 



0.04 ±0.01. 



(22) 
(23) 



The calculated values of A'o and fo both clearly have a linear 
dependence on the assumed value of EMo . The estimates for 
A^o and fo will correspondingl y change if a slig htly different 
value of EMo is assumed (e.g., lHill et al.ll2008f) . but the val- 
ues for all other quantities discussed above (i.e., no, DMo, 
Hn, Hn2, Hn & Hf) will be unchanged. 

For the sake of completeness and for comparison with 
previous studies, we also compute the vertical "occupation 
length", Lc = (DMo)^/EMo, and "characteristic de nsity", 
Nc = EMo /DMo, of the WIM (Reynolds 1991b; Hill et al.l 
I2OO8). For DMo ~ 26.0 pc cm"^ and EMo « 1-9 pc cm"® 
as derived above, we find Lc « 350 pc and Nc « 0.07 cm"^. 
However, we caution that in the case where Hn/Hn2 7^ 2 as 
argued here, these quantities have no direct physical interpre- 
tation — see further discussion in i|4.4l 



4 Discussion 

Our main result is that within a few kpc of the Sun, the scale 
height of the thick disk of warm ionised gas in the Galaxy is 



'"Note that the fitted value of no is not the true mid-plane density, 
but that of the thick-disk component alone, extrapolated to z = 0; 



see discussion in i]4.4l 
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H„ ~ 1.8 kpc. This is a factor of ~ 2 higher than the value 
that has been found in most previous studies (see summary 
m Table [B- If this is a global property of the Galactic ISM, 
this has a number of important implications. 

However, it is important to note that our large value for 
H„ results from our exclusion of low-latitude pulsars, whose 
DMs require a higher value of no, and hence a lower value 
of H„, for constant DMq. Because our sample consists only 
of high-|6| pulsars, it could be argued that we are probing the 
WIM immediately above and below the Sun, and thus may 
be tracing local stru cture rather than an y ov erall properties of 
the Milky Way (see lRevnoldsHl97lfl991ah . 

In the following discussion, we first consider the possible 
presence of local structures in DM and in EM, and conclude 
that our results are not biased by any such features. Next, 
we examine our implicit underlying assumption that a single 
phase of the WIM simultaneously produces both DMs and 
EMs, and confirm that ths assumption is justified by the data. 
In the rest of this section, we then consider some more gen- 
eral implications for the properties of warm ionised gas at 
increasing distances from the Galactic mid-plane. 

4.1 The Effects of Local Structure 

4.1.1 Local Structure in Dispersion Measure 

In Figure|4]we show the locations of the pulsars in our sample 
projected onto the Galactic plane, overlaid with the NE2001 
electron density model. Of the 15 systems for which |fe| > 
40° (shown in blue), six sight lines (all to globular clus- 
ters) extend over a significant distance through the Galaxy, 
traversing several, separate, spiral arms. These data clearly 
do not probe merely local gas, but are representative of the 
WIM within several kpc of the Sun. In any case, these sight- 
lines should not be dominated by individual density features, 
since the long path lengths ensure that individual clumps or 
cavities in the electron distribution largely average out. 

The remaining nine systems (all pulsar parallaxes) lie 
much closer to the Sun. As indicated in Figure ID there 
are severa l cavities and enhancements in the solar neighbour- 
hood (see ICox & Revnoldsll987l : lToscano et al.ll999l : ISfeir et al.l 
1 19991 . and references therein), whose presence could domi- 
nate the DM of nearby pulsars, and thus might bias our fits. 
A qualitative argument against this possibility is that in Fig- 
ure[T] the envelope defined by the high-|6| parallax data joins 
smoothly onto that traced out by the high-|6| globular cluster 
data. Since we have argued above that the DMs to globular 
clusters trace the widespread properties of the WIM in the 
disk, the pulsar parallaxes are also consistent with this inter- 
pretation. 

We can quantify the effects that local structure might be 
having on the data by plotting DM/d versus z, where d is 
the distance to each source. This distribution is plotted in 
Figure [5] the electron column per unit distance plotted in 
the ordinate is equal to the mean electron density along the 
sightline to a given source. The corresponding predictions 
for DM(2)/d = DMx(2) /z for each of the curves drawn in 
Figure [T] are plotted in Figure [5] with the form of DM±{z) 
given by Equation ((SJ. 

For a generally smooth distribution of electrons in all di- 
rections, values of DM/d will follow a single curve such as 
those plotted in Figure [5] Pulsars located inside or just be- 
yond a cavity will exhibit values of DM/d that sit below this 



curve, while pulsars that sit behind density enhancements 
will have values of DM/d above this curve. Both such ef- 
fects are clearly seen in Figure [S] especially for sources at 
low \b\ (coloured red and orange), which are not included in 
our best fit. However, if local structure is affecting our data, 
then high-latitude pulsars (coloured blue) at low values of z 
should have values of DM/d that deviate more from the best 
fit curve than do data for pulsars at higher z. Only for one 
pulsar (at z ~ 300 pc) of the 15 sources in our sample is 
there a slight suggestion of such an effect. The small errors 
in no and Hn derived from bootstrapping confirm that these 
data do not have a significant effect on our overall fit. 

From Figures |4] and |5] we conclude that localised struc- 
tures are not biasing our results, and that the fits to DM data 
in Figure [T] are a meaningful representation of the ensemble 
scale height of diffuse ionised gas within several kpc of the 
Sun. 

4.1.2 Local Structure in Emission Measure 

Another way in which localised structure may bias our mea- 
surements is in our assumption that the scale heights //„2 
and Hn co rrespond to the same regions. The determination 
of Hn2 by H affner et alj ( Il999h was specific to a particular 
segment of the Perseus arm. In contrast, the DM data in Fig- 
ure[T]represent integrated column densities spread over a fac- 
tor of ~ 200 in z and distributed over a range of Galactic 
longitudes (see Figs.|2]and|4j. Thus while we have argued in 
ij4.1.1l that the DM data represent ensemble properties over 
several kpc of the Galactic disk, the EM data provide no such 
guarantees. If ionised gas above the Perseus arm has special 
properties (and indeed Fig. |4] shows that none of our high- 
latitude pulsars lie in or behind the Perseus arm), then the 
conclusions of ii3.3l are not correct. 

There are two arguments that support the calculations in 
il3.3l First, Lazio & Cordes (1998a,b) have specifically stud- 
ied the properties of diffuse ionised gas toward the Perseus 
arm and over the rest of the outer Galaxy, and conclude that 
the ionised disk shows neither warping nor flaring in these 
regions. Second, even if structure peculiar to the Perseus 
arm might introduce some bias into our calculations, the fun- 
damental result in Sj33]is that //„/J/„2 2> 2. Only com- 
pletely different behaviour for the Perseus arm compared to 
the sightlines in other directions probed by pulsar DMs could 
reduce this ratio sufficiently to remove the need for a z-dependent 
filling factor. NGC 891 and M31, two spiral galaxies with 
broadly similar properties similar to the Milky Way, show 
variations in the thickness of their ionised gas layers as a 
function of position, but at fract ions of no more than ~50% 
tRand et al.. 1990 : Miller & Veilleux 2003: Fletcher et al.l2004l) . 

4.2 Two Cloud Populations? 

In ii3.3l we assumed that at a given height z, WIM clouds 
can be described by a single volume-averaged electron den- 
sity, n{z). This is obviously an idealised situation, but is a 
reasonable approximation of actual data, which suggest that 
n{z) shows a modest spread in density around a well-defined 
average value (ffill et aljf2008b . In the model considered in 
ii3.3l we have assumed that this simple cloud population pro- 
duces both the DMs and EMs that we observe. The fact 
that Hn/Hn2 7^ 2 can then be explained by invoking sep- 
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arate dependences of A'^ and / on z, as described by Equa- 
tions Gil & i fTB . 

However, an alternative possibility is that there are (at 
least) two populations of WIM clouds, with separate densi- 
ties, filling factors and scale heights. Specifically, if we sup- 
pose that one set of clouds has a low density and a high filling 
factor, while a second cloud population have a much higher 
density and a much lower filling factor, then observations of 
DM will only probe t he first group , while data o n EM will 
only trace the second ('Gould l 197ll : iHeileslllOOll) . Carrying 
out a joint analysis of DMs and EMs to infer overall proper- 
ties of A'^ and /, as we have done in ij4.2l would then not be 
meaningful. 

We thus now consider whether the data can be validly in- 
terpreted in terms of two separate populations of WIM clouds. 
Suppose that one group of WIM clouds has electron den- 
sity Na{z), filling factor /a (which is not a function of z) 
and volume- averaged density nA{z) = fANA{z). The den- 
sity of these clouds is distributed exponentially, such that 
nA{z) and Na{z) both have a scale height Ha- A second 
cloud population has corresponding parameters Nb{z), /s 
and 713(2) = fENsiz), with an exponential scale height 
Hb. 

The vertical dispersion and emission measures integrated 
through the thick disk then become: 

DMo = fANA,oHA + fsNBfiHB , (24) 
EMo = fANloHA/2 + /s7V|,oHs/2 , (25) 

respectively, where the subscript "0" indicates extrapolated 
mid-plane values, and where we have disregarded second or- 
der terms in /a/s. The observed behaviours of DM and 
EM with 2 are both well described by a single exponential 
layer (Fig. \T\ Fig. 4 of Haffner et al. 1999), implying that 
only a single population of clouds contributes to each set of 
measurements. We thus adopt the requirement that the DMs 
only probe population A but that the EMs only probe popu- 
lation B: 



4.3 General Implications 



and 



fBNBfiHB < DMo 



fBNl^oHB/2 « EMo. 



(26) 



(27) 



Comparing Equations Olt & l l27t . we can then identify Hb /2 
H^2, so that combining Equations l l26t & J27t then vieldsPI 

L. 



fB « 



iH„ 



0.2, 



Nb,o > 2Nc ~ 0.15 cm-^ 



(28) 
(29) 



with Lc and Nc defined at the end of §33] We can rule out 
the existence of such a component to the diffuse WIM, be- 
cause Hq data lead to direct estimates of mid-plan e densities 
for the diffus e WIM of Nb o « 0.1 - 0.3 cm~^ dRevnoldsl 
[l980. 1983; Revnolds et al.ll995l) . Since the predicted cloud 
densities of population B are at odds with those observed, we 
conclude that it is not valid to argue that DMs and EMs are 
produced by separate populations. The alternative is that the 
distribution of ionised clouds is described by a single set of 
parameters, 71(2) = f{z)N{z), as laid out in §33] 



"a similar analysis for population A implies fA ^ 0.1 and 
Nao <■ 0.15 cm-3. 



We have shown in ^4. 1 1 that our results are not likely to be 
merely a result of local structure in DM or EM, and we have 
demonstrated in §42] that we have not mistakenly conflated 
data on two separate populations of ionised clouds. We em- 
phasise that even if some of our arguments in these Sections 
can be shown to have shortcomings, our main result, that 
Hn « 1.8 kpc, is derived purely from pulsar DMs, and does 
not rely on any subsequent assumptions about filling factors 
or EMs. We now consider some of the implications of our 
revised WIM parameters. 

We first consider the Galactic distribution of pulsars. The 
most frequent application of models for the WIM is to de- 
rive distance estimates of pulsars from measurements of their 
DMs. Figure [T] demonstrates that the NE2001 model does 
not correctly reproduce the DMs of high-latitude pulsars at 
known distances. In turn, it seems likely that the DM-derived 
distances for many high-|6| pulsars without independent dis- 
tance constraints may be substantial underestimates; the cor- 
responding pre dictions from NE200 1 s hould be used with 
caution. Indeed. iKramer et^ ( l2003l) and lLorimer et aLl ( l2006h 
found that if distances predicted by NE200I are applied to 
the Parkes multibeam pulsars, these sources are then anoma- 
lously clustered toward small 2. 

A revision of the distances to high-latitude pulsars would 
have important implications. Such data play a major role in 
characterising the population statistics of millisecond (which 
tend to be found at high \b\ due to their large ages), as well 
as to predict the yields of future pulsar surveys. Recent stud- 
ies of these i ssues have all been underpinned by th e NE2001 
model (e.g.. iHobbs et al.ll2005l; Istorv et alj|2007l) . and the 
corresponding conclusions may now need to be reconsidered. 
Specifically, if the true distances to high-|6| pulsars are larger 
than have been assumed in recent studies, then millisecond 
pulsars may have higher space velocities, be more luminous, 
and be more numerous than has been previously argued. 

Overall, we note that while the NE2001 electron den- 
sity model most likely underestimates the distances to high- 
latitude pulsars, one of the major strengths of this model was 
that it corrected the severe over-estimates that previous mod- 
els had made for the distances to /ow-latitude pulsars. Thus it 
is unwarranted to argue that the entire NE2001 model is un- 
reliable; our analysis focuses only on the high-latitude com- 
ponent, and indicates only that this aspect needs revisiting. 
Correspondingly, we do not recommend that Equation l|6j or 
Figure [T] be used to estimate the distances or DMs to indi- 
vidual pulsars at high |fe|, but instead suggest that in future 
studies, the results derived here should be incorporated into 
a revised global model for the disk, spiral arms and halo. 

As noted in §TTTTT] a full three-dimensional model of 
Galactic magnetic fields, cosmic rays and thermal electrons, 
but assuming a warm layer with scale height H„ ~ 1 kpc, re- 
sults in an unacceptably high halo magnetic field strength of 
~ 10 fiG, and an u nreasonable trun cati on of cosmic ray elec- 
trons at 2 = 1 kpc dSun et al.l200§) . As lSun et al] ( l2008» ex- 
plicitly discuss, a WIM layer with a scale height of ~ 2 kpc 
can simultaneously solve both these problems. The indepen- 
dent conclusions that we have arrived at here thus support a 
relatively weak (~ 2 ^G) magnetic field in the halo, and a 
cosmic ray distribution which can diffuse to high 2. 

Despite the more consistent overall picture of the WIM 
that our analysis can potentially provide, we note that even 
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after artificially increasing the errors in z as was carried out 
in i]3.1.1l the best fit in Figure [3] has ~ 5. The fact that 
Xr S> 1 clearly indicates that a smooth exponential slab 
is not an accurate description of the individual data points. 
As discussed extensively by Cordes & Lazio (2002, 2003), 
many sightlines through the Galaxies pass through interven- 
ing regions where the density is either enhanced (clumps) 
or reduced (voids). These systematic deviations are not in- 
corporated in our model, and imply that Xr will always be 
relatively large in any simple fit to this complex data set. 
However, it is noteworthy that Xr does drop significantly for 
higher cut-offs in |6| (see Fig. [3j, indicating that the WIM 
becomes more homogeneous as one moves further from the 
Galactic plane. We further note in Figure [T] that the low- 
latitude data (in red) are distributed both above and below 
the fits that include these data (dashed and dot-dashed lines), 
while the highest latitude data (in blue) generally lie only on 
or above the corresponding best fit (solid line). This suggests 
that in the Galactic plane there are both clumps and voids in 
the WIM, but that far from the plane any deviations from a 
smooth distribution are mainly in the form of clumps, rather 
than voids. This same argument produces additional quali- 
tative evidence against the NE2001 model for ionised gas at 
high z. If the NE2001 description (shown as a dot-dashed 
line in Fig.[T} is correct, then the fact that almost every blue 
data-point in Figure [T] lies well below this line requires that 
there be separate voids in electron density in front of all these 
high-|6| pulsars. A large number of voids to account for all 
these low DMs does not seem warranted, especially given 
that these pulsars are distributed at a wide range of distances, 
in a variety of directions, and sit both above and below the 
Galactic plane (see Figs. |2] & |4j. The sustained low DMs 
for all these pulsars compared to the NE2001 model suggests 
that a different overall density distribution is required in these 
directions. 

4.4 Mid-Plane Properties of the Diffuse 
WIM 

In this section, we compare our extrapolated mid-plane val- 
ues for the thick-disk WIM to other studies of the ISM in 
the disk. At 2 = 0, we infer a volume-averaged density for 
the diffuse WIM of no = 0.014 ± 0.001 cm"^. Observa- 
tionally, this can be compared to values of DM/d for pul- 
sars at low |6| with known distances, which typically yield 
no ~ 0.03 — 0.1 cm^^ along complicated sightlines through 
spiral arms, and np ~ 0.01 — 0.02 cm^*^ in i nter-arm regions 
dJohnston et 311 120011 : IWeisberg et alj|2008l and references 
therein). As discussed in !jl] this indicates two components 
for warm ionised gas in the ISM: free electrons associated 
with individual H II regions (only found in spiral arms at low 
|fe|), and the diffuse, thick-disk WIM (seen at all latitudes 
and in all directions). We have not included any contribu- 
tion from the former component in our model, since we only 
fit DMs at high \b\ and high z. Thus the fact that observa- 
tional estimates for no = DM/d in inter-arm regions agree 
well with our derived value for no supports the picture of an 
extended diffuse WIM in which individual dense ionised re- 
gions are embedded, and suggests that the thick WIM layer 
extends all the way down into the Galactic plane. 

Estimates of no from previous attempts to globally model 
the thick-disk WIM are summarised in Table [T] Our value 
for no is systematically lower than that in almost all these 



other studies. In most cases, this results from the fact that 
DMo = noHn for an exponential slab. Since all studies 
generally agree on DMo, but most previous authors have de- 
rived much smaller scale heights than we have found, their 
estimates of no are correspondingly larger. BMM06 set only 
a lower limit on //„: their estimate of no may be biased by 
the systematic distance uncertainties in their analysis, as dis- 
cussed in i] 1.1. 21 

A variety of studies have also estimated the mid-plane 
filling factor of the WIM, as summarised in Table |2] The 
value fo — 0.04 ±0.01 which we derived in i]3.3l agrees with 
the estimate of BMM06, but is substantially smaller than that 
found in most other studies. This discrepancy results from 
the fact that calculations of fo are us ually degenerate with 
estimates of no (see, e.g., Eqn. [9] of iRevnoldsl 1 1 99 1 bh . In 
our model, fo oc (no)^ (from Eqn. II23I ). and we have ex- 
plained above that no needs to be revised downward from 
previous estimates by a factor of --^ 1.7 — 2. The true value 
of fo in the WIM must then be a factor of 3 — 4 smaller 
than what has been found in previous studies, as we have 
determined here. 

Our estimates for fo and no imply a typical mid-plane 
density for WIM clouds iVo = no /fo = 0.34 ± 0.06 cm"^, 
w hich agrees with the theoretical prediction No = 0.3 cm~^ 
of lCoxl ( |2()05|). Independent estim ate s of Np have been car - 
ried out bv lRevnoldsl l l 1 980LI 1 983h and lRevnolds et al.l ( Il995h . 
who measured EMs along a variety of sightlines and found 
A^o ~ 0.1 — 0.3 cm^'', in reasonable agreement with our 
value. The thermal pressure of the WIM iQ f = 2NkTe, 
so for Te ~ 8000 K (e.g., Reynolds et al. 1999), we infer a 
mid-plane ISM pressure for the thick-disk WIM of Po/k ~ 
5400 ± 1000 K cm~^. Other estimates for the pressure of 
the WIM and of other ISM phases show a wide variation de- 
pending on local conditions, but certainly the above value 
falls comf ortably within the typically observed range (see 
ICoxll2005l) 

iHeilesI ( 1200 ll) argued that there must be two separate 
components of the diffuse WIM, one which contributes to 
EMs seen in Ha, and a second which separately produces the 
DMs seen for pulsars. We have already shown in i]4.2l that. 
as an overall model for the ISM, a two-phase WIM is not 
supported by the da ta. Here we br iefly consider the specific 
arguments made by iHeilesI ( 1200 ih . He begins by assuming 
that the EMs and DMs are produced in the same ISM phase, 
and then writing No = Nc = EMq /DMq. The implied den- 
sity leads to an anomalously underpressured WIM compared 
to the cold, neutral, ISM, a problem which can be relieved 
by putting the electrons that contribute to EMs and DMs into 
distinct phases. However, the expression A'^o = EMo /DMo 
is only valid for a slab of uniform density; for the expo- 
nential fall-off considered here, the appropriate expression 
is iVo = (EMo/-ff„2)/(DMo/-ff„), which leads to a WIM 
pressure ~ 4 — 5 times higher than in the discussion of lHeilesI 
(2001). With this correction, there is approximate pressure 
equilibrium between the WIM and other phases of the ISM, 
so that a model in which the EMs and DMs are produced by 
the same gas is sustainable. 

Finally, we can combine Equations ( I16t & ( I21t to esti- 
mate a mean square electron density at midplane (n^)o = 
EMo///„2 = 0.005 ± 0.001 cm~^. This value, estimated 
from quantities which we have not substantially revised in 



The factor of 2 coiresponds to a fully ionised hydrogen gas. 
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this paper, is broadly comparable to previous estimates, as 
listed in Table[3] 

In summary, our results imply that within a few kpc of 
the Sun, the thick-disk component of the WIM at mid-plane 
is comprised of clouds of electron density A'^o ~ 0.3 cm^'^ 
and a volume filling factor of /o ~ 4%, corresponding to 
a volume-averaged electron density no ~ 0.014 cm~'^. If 
the diffuse WIM is considered separately from denser gas 
in the plane associated with discrete H II regions, these val- 
ues are in broad agreement with independent estimates made 
by other approaches. We note that radio recombination line 
studies generally identify higher values of A'o a nd lower val- 
ues of fo than what we have derived he re fe.g.. lHeiles et al.l 
ll996bl ; lRoshi & Anantharamaiahll200lh . but these emission 
lines preferentially trace denser gas, possibly associated with 
the outer envelo pes of H II regions or w i th the walls of chim - 
neys in the ISM jAnantharamaiaUl986l : lHeiles et aljl996bh . 



4.5 Anti-Correlation of and /, and Evo- 
lution of WIM Properties with z 

The total mass and ionisation rate of the WIM are u nchanged 
in o ur analysis compared to previous estimates fe.g.. lRevnoldsl 
[T990b), since the calculations of these quantities primarily 
depend on the integrated column densities EMo and DMo, 
rather than on the scale-height. What the calculations in 
fj33]do imply is contrasting behaviours for the cloud den- 
sity, A'', and the locally averaged electron density, n, as a 
function of height above the Galactic plane. The former de- 
cays over a relatively small scale height, ~ 500 pc. But the 
volume filling factor, /, of these ionised clouds increases as 
one moves further from the plane, counter-balancing the de- 
creasing internal cloud density. Since / increases with z 
at a rate slightly more slowly than that with which A'^ de- 
cays, n = fN decays relatively slowly as a function of 
z. These results agree with the earlier findings of BMM06, 
who through a separate approach derived broadly compara- 
ble values of Hm and Hf to those found here, albeit with 
much larger uncertainties (compare Eqn. O with Eqns. 1191 
&|20l). 

An important distinction between our analysis and that 
of BMM06 is that we find Hn < Hf, while BMM06 con- 
cluded that Hn ~ Hf. As per Equation {Hll, BMM06 thus 
could infer only a lower limit on H„, and indeed they argued 
that there is no evolution in n for < z < 1 kpc, based 
on the constancy of DM/d over this range. However, as dis- 
cussed in i ]1.1.2l there are large uncertainties in the values 
of d and z used in their analysis. Figure [5] shows that for 
our sub-sample of pulsars with reliable distances and with 
|6| > 40°, there is a gradual decline in DM/d with increas- 
ing z up to 2 ~ 1 kpc, which joins smoothly to the behaviour 
seen for z ^ 1 kpc (note that the scatter for the blue points in 
Fig.|5]is much smaller than for the equivalent plot in Fig. 12 
of BMM06). This can only occur if n steadily decays with 
z, requiring Hm < /// as we have inferred above. 

The decaying exponential for A'^ compared to the grow- 
ing exponential in / obviously implies an overall anti-correlation 
between / and A' ^, as has also been found through various 
other approaches ( iPvnzaril993l . lBerkhuiisenll998l : lBerkhuiisen < 
|2008| . BMM06). For a dependence: 



(30) 



we can combine Equations ([To), Ol, GS), (|20), ^ & l|23ll 
to derive a = -0.73l[J;|5^ and (3 = 0.018 ± 0.004. This 
determination differs substantially from a — —1.07 ± 0.03 
found by BMM06, but agrees with a ~ —0.7 as derived by 
iPvnzar'l (1993). Since a = —HN/Hf, the discrepancy be- 
tween our estimate and that of BMM06 simply results from 
our finding that Hn < Hf as discussed above. 

BMM06 consider two possible p hysical interpretations 
of Eq uation ( 130b (see also discussion bv lBerkhuiisen & MuUerl 
120081) : either it can represent a stochastic, turbulent ensem- 
ble of clouds of different sizes and different densities, or it 
describes the equation of state of ionised clouds, in that it 
probes the response of the WIM to a change in external con- 
ditions. In the situation that we consider here, the variation 
of / and A'^ are both modeled as smooth functions of z, so 
the latter case applies. 

In this context, what is the physical significance of the 
anti-correlation between A'^ and /, and of the differing scale 
heights for these two quantities? We first note that the ther- 
mal gas pressure, P = 2NkTe, depends o n A' but not on 
/. If the evolution of Tg with z is modest jRevnolds et al.l 
119991 : iPeterson & Webbed 1 20021) , then the pressure decays 
exponentially with a scale height Hn ~ 500 pc, as de- 
scribed by Equation dl It . This agrees well with the pres- 
sure scale height derived from hydros tatic equilibrium calcu- 
lations of the f ull multi-phase ISM feoulares & Co}3ll990l : 
lFerrierell998ah . 

To interpret the evolution of / with z, we consider the 
simple case in which the WIM is composed of a large number 
of clouds, all with the same number of electrons per cloud, 
A/", and with Cdz clouds contained in a vertical interval dz. 
The volume of an individual cloud is V{z), so that A'^(2) = 
M /V[z). If the total surface area of the Galaxy (assumed 
to be independent of z) is E, then /(«) = CV{z)/Tj — 
{CN /T,)N{z)~'^ . The exponential decay of Af(2) must then 
imply matching exponential growth in V{z) and in f(z): 
clouds will be proportionally larger in regions where the am- 
bient pressure is lower. This case corresponds to a = —1 
(i.e., Hn = Hf), as claimed for the WIM by BMM06. 

We can now generalise the above scenario, and allow C 
and M to both evolve with z. If q = —1, then J\f{z) oc 
C[z)^^ , i.e., individual clouds atz = z\ might contain more 
electrons than those at 2 = 22, but in this case the total num- 
ber of clouds at 21 will be smaller than that at 22 by the same 
factor. The total mass of the WIM remains constant as a 
function of 2. 

Now consider the case we have derived here, with a = 
-0.73 (and Hn < Hf). The fact that a > -1 indicates that 
a cloud of a given electron content will be larger in a region 
of low pressure compared to one under high pressure, but 
that the increase in size will be less than that expected from 
the ratio of pressures alone. This could be due to Te decreas- 
ing with 2, but this would, if anything, be opposite to the 
observed trend (Revnolds et al. 1999). The alternative is that 
C[z)M{z) is decreasing with 2, i.e., the total particle content 
of the ensemble of clouds (and hence the mass of the WIM) 
at high 2 is smaller than for those at low 2. Quantitatively, we 
can write n{z) = f{z)N{z) = C{z)M{z)/'L. Equation ^ 
T^c^^iiiplies that N {z)C {z) oc e"^''^". We interpret this as 
n;ne mcreasing scarcity of warm gas at higher 2, tracing the 
transition from the WIM into the hot halo. Electrons at the 
surface of WIM clouds are heated to 10^- 10"^ K via thermal 
conduction jCowie & McKed 19771 : 1 Vieser & Henslej2007h ; 
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they then no longer emit in Ha, and have too low a density 
to be seen in pulsar dispersions (see i]4.6l below). Such ma- 
terial is thus no longer observable through the DMs or EMs 
considered in ijs] 

The implication of the above discussion is that of the var- 
ious parameters (Hn, H„2 , Hn, Hf, a, f3) that we have used 
to describe the behaviour of the WIM as a function of z, the 
two quantities that best encapsulate the physical conditions 
of the gas are Hn and Hn. The former traces the approxi- 
mate fall in gas pressure away from mid-plane, as moderated 
by the gravity of the disk and by the effects of magnetic and 
cosmic-ray pressure. Separately, the latter reveals the extent 
to which the WIM permeates into the halo, and so probes the 
lifetimes of warm clouds embedded in hot coronal gas. Since 
these are two separate physical processes, it is unsurprising 
that they have distinct scale heights. The fact that Hn 7^ H„ 
then demands the existence of a third scale height, Hf, which 
describes the geometry and size that clouds adopt as they ac- 
commodate external pressure but at the same time evaporate. 

4.6 Constraints on Hot Coronal Gas 

At significant distances from the Galactic plane (z 2000 pc), 
theory predicts that Z 90% of the ISM b y volume should 
be occupied by lO'^ - 10*^ K coro nal gas ( lFerrierel[T998bl: 
iKorpi et al.|[l999l : Ide Avillej|200(]h . To interpret our data 
at these heights above the plane, we first must determine 
whether free electrons associated with c oronal gas can con- 
tribute significantly to the observ ed DMs jWolfire et aljl995l ; 
iBregman & LIovd-Daviesl2007t) , or whether we are only prob- 
ing the WIM. 

The standard picture of the Galactic halo is an extended 
high- 2: region of v ery low electron density, nhaio ^ 10~^ — 
10"* cm"^ (e.e.. !Sembach et alj|2003l : I Yao & Wane 2005h. 
We can test for the presence of such gas in our DM data by 
attempting a linear fit of DM. sin |6| vs. z to the data in Fig- 
ure [T] (and relaxing the previous constraint that |6| > 40°), 
but only considering systems for which z > 5000 pc, where 
the contribution from the WIM should be small. The best 
fit slope to these six measurements is consistent with zero; 
bootstrapping provides a Sa upper limit on the halo density 
nhaio < 7.6 X 10^* cm^''. We conclude that any contribu- 
tion to pulsar DMs from hot, low density gas in the halo is 
below our detection limit. 

An independent insight into the relative contributions of 
the WIM and corona to pulsar DMs comes from the study of 
Howk et al. (2003, 2006), who compared a variety of ISM 
tracers, including pulsar DMs, in the foreground of the glob- 
ular cluster NGC 5272 (Messier 3). NGC 5272 contains 
four radio pulsars, and corresponds to the blue data point 
with z = 10.2 kp c, DM. sin|fe| = 25.9 pc cm"^ in Fig- 
ure[T] iHowk et al.l ( 2006) compute separate electron column 
densities along this sight line for the WIM, for ~ 10"" K 
coronal gas and for ~ 10® K coronal gas. The total DM 
inferred from pulsa r observati ons is 26.4 ± 0.3 pc cm~^ 
dHessels et aDl2007! ). of which iHowk et all j2006h estimate 
that the WIM contributes ~ 25 pc cm~ while 10^ K gas 
contributes ^ 1.5 pc cm"''. Any remainin g contribution 
from 10® K gas must be small: iHowk et aD (2006) adopt a 
model in which this gas has an average electron density of 
2.4 X 10"'' cm"^. The DM contribution from the globular 
cluster itself is ex pected to be negligible dFreire et alj|200ll : 
iHowk et alj|2003l) . These results confirm that almost the en- 



tire electron column out to high z originates in warm gas, 
at an electron temperature ~ 10* K. The pulsar DMs in Fig- 
ure[T|thus predominantly trace the WIM, and do not constrain 
the properties of hotter, coronal gas. 

4.7 Warm Ionised Gas in the Galactic Halo 

Given that our data do not appear to be substantially contam- 
inated by coronal gas, we can now consider what our results 
imply for the filling factor and pressure of the WIM at high z. 

The physical requirement that / < 1 implies that Equa- 
tion dlOt and the results that follow are not well-defined for 
heights z > 2200 ± 300 pc. The value i?„2 = 400 ± 30 pc 
has been derived from Ha data in the range 6° < |fe| < 35° 
jHaffner et al .11 19991) . which at the distance to the Perseus 
arm corresponds to a maximum height z = 1400 pc. The 
EM data used to derive the dependence of / on z thus fall in 
the regime where / < 1, so that our approach is internally 
self-consistent. At the top of this layer, we find / = 0.30, 
TV fa 0.022 cm"^ and n = fN ^ 0.0065 cm'^ This 
agrees with the findings of BMM06, who argue that f{z) 
reaches a maximum value / « 0.3 at 2; « 1 kpc. How- 
ever, this calculation contrasts with previous claims that the 
WIM is the dominant p hase of the ISM at 2 ~ 1 kpc (e.g., 
iRevnoldsl [19891 Il991ah . While the WIM probably domi- 
nates over neutral gas in these regions, at this height coronal 
gas already occu pies most of the volume of the ISM (e.g., 
.Korpietal .1119991) . 

Beyond the top of the WIM layer studied by WHAM in 
the Perseus arm. Equation i ll It cannot be meaningfully in- 
terpreted, since we no longer have independent constraints 
on both n and n^. However, we can broadly characterise the 
pressure and filling factor of the WIM at high z by noting 
that P oc N oc n/ f. At large z, we know that n continues 
to decay at least as fast as the exponential in Equation (|5]l, 
because of the high quality of the fits of this model to the 
data in Figures [T]& [5] (for z ^ 5 kpc, we cannot rule out the 
possibility that n begins to fall off more rapidly). While / 
appears to increase with distance from the mid-plane at low 
values of z, this function must turn over at 2 2000 pc to 
accommodate coronal gas (e.g., Korpi et al. 1999l: lde Avine3 
I200Q) . If n and / are both decreasing functions at high z, 
then A'^ = n// at these heights must become increasingly 
larger than predicted by Equation i ll lb . This then implies 
that any warm ionised gas at high 2 is over-pressured relative 
to the exponential decay in pressure seen within 2 kpc of the 
mid-plane. This high pressure is consistent with the expec- 
tation that any small warm c louds at these he ights should be 
enveloped in a hot halo (e.g.. lFox et al ]l2005h . 

At 2 ~ 5 — 10 kpc, the ambient pressu re in the Galac- 
tic halo is P/k fa 1000 - 2000 K cm"^ I'S huU & Slavinl 
1994; Sternberg et al. 2002; Everett et al. 2008). If we as- 
sume 10* K for the WIM, and adopt Equation 1(5) as an 
upper limit on n at these values of 2, we infer / < 0.02 for 
any warm ionised gas at these heights above the plane. These 
small cloudlets may correspond to th e "drizzle" of ionised 
gas a ccreting into the halo proposed bv lBland-Hawthorn et alj 
(2007), or to the extraplanar droplets and clump s of Ha seen 
in observations of some edge-on spiral galaxies l iFerguson et alj 
ll996l ; ICeciletalj200lh . 
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5 Conclusion 

We have used the dispersion measures of pulsars at known, 
reliable, distances and the emission measures traced by dif- 
fuse Ha emission to re-evaluate the structure and distribution 
of the thick disk of warm ionised gas in the Milky Way. Our 
main conclusions are as follows: 

1. The DMs of high-latitude pulsars are well fit by a 
plane-parallel distribution of free electrons, in which 
the volume-averaged electron density, n, decays ex- 
ponentially with height, z, above the Galactic plane. 

2. Within a few kpc of the Sun, the exponential scale 
height for n is Hn ~ 1800 pc, approximately double 
the value determined in most previous studies. Esti- 
mates of distances to high-latitude pulsars using the 
"NE2001" electron model of Cordes & Lazio (2002, 
2003) are thus likely to be substantial underestimates. 

3. Within a few kpc of the Sun, the mid-plane volume- 
averaged electron density of the diffuse WIM is no — 
0.014 ± 0.001 cm~^, and the corresponding total ver- 
tical column density of the WIM is ~ 26 pc cm~^. 
Assuming a total vertical emission measure EMo — 
1.9 ± 0.3 pc cm^'', clouds in the thick-disk compo- 
nent of the WIM have a typical internal electron den- 
sity of iVo = 0.34 ± 0.06 cm"^ at mid-plane, with a 
volume filling factor fo = 0.04 ± 0.01. 

4. Comparison of DMs and EMs shows that the inter- 
nal electron density, A'^, of individual ionised clouds 
within the thick-disk WIM decays rapidly with z, with 
an exponential scale height Hn ~ 500 pc. The dif- 
fering behaviour for Hn and Hm requires that the 
filling factor of the WIM grows exponentially with z 
over the range < 2 ^ 1.4 kpc, with a scale height 
Hf ^ 700 pc. 

5 . The filling factor of the thick-disk WIM probably reaches 
a maximum of / ~ 0.3 at z ~ 1.5 kpc. At larger 
distances from the Galactic plane, the ISM becomes 
increasingly dominated by hotter, coronal, gas in the 
halo. However, hot low-density halo gas makes only 

a negligible contribution to observed pulsar DMs. 

A variety of future observations will allow us to test the 
above conclusions, and to better understand the structure of 
the diffuse WIM. Most immediately, from late 2008 WHAM 
will be relocated to the southern hemisphere, allowing stud- 
ies similar to those carried out on the Perseus spiral arm to 
now be repeated for other regions. 

Further improvements will come from an increase in the 
number of high-latitude radio pulsars with direct distance es- 
timates. For example, only ~ 25% of the Galactic globular 
clusters at |6| > 40° contain known radio pulsars. While 
some clusters at these latitudes have been searched for pulsa- 
tions without succe ss (Anderson, 1993; Biggs & Lvne.1996 ; 
iHessels et aLll2007h . at least some of these target s probably 
conta in as yet undetected pulsars (e.g., NGC 5634: ISun et al.l 
l2002h . and should be searched more deeply. Furthermore, 
there are several high-latitude pulsars with 1.4 GHz radio 
fluxes at the level of a few mjy, which is sufficiently bright 
to derive interferometric parallaxes through future observa- 
tions. 

Most surveys for radio pulsars have concentrated on the 
Galactic plane, and those surveys that have covered high- 



latitude regions have been of co mparatively modest sensitiv- 
ity (e.g., [Manchester et aljl99^ . Forthcoming all-sky pulsar 
surveys with more sensitive instrumentation will greatly ex- 
pand the sample of high-|6| pulsars, whose DMs can then be 
used to better map out the structure of the WIM. 

Finally, we note that the new model we have presented 
for the WIM can also be used to probe Galactic magnetic 
fields at high z. In a subsequent paper, we will combine the 
density distribution derived here with Faraday rotation mea- 
sures of a large number of polarised extragalactic sources at 
high with which we aim to derive the strength and geom- 
etry of magnetic fields in the Galactic halo. 
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Figure 1 : Distribution of dispersion measure as a function of height above the Galactic plane for pulsars with 
known, reliable, distances. Pulsars with distance determinations from trigonometric parallaxes, from associations 
with globular clusters and from associations with the Magellanic Clouds are shown as filled circles, open circles and 
stars, respectively (data for the Magellanic Clouds are indicated as upper limits, since the Galactic and Magellanic 
DM contributions cannot be separated). Data are colour-coded according to the magnitude of their Galactic latitude, 
|6| , as indicated in the legend. The error bars represent the published uncertainties. The three curves represent planar 
slabs of ionised gas with a scale height i7„ for the volume averaged electron dens ity, n. T he dot-dashed line shows 
a n cx se ch^ (z/Hn) d ensity dependence, using the parameters given in Table 3 of lCordes & Lazio (2002) (see also 
Fig. 1 of ICordes & La zio 2003); the dashed line shows the best fit for an exponential slab (n oc e~^/^") to all the 
data in the figure; the solid line shows the best exponential fit to only those data for which \b\ > 40°. 
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Figure 2: Aitoff projection in Galactic coordinates of pulsars with know n, reliable, distances, super-imposed on 
the all-sky distribution of Ha intensity (as derived by lFinkbeinerll200"3h . The Galactic Centre is at the centre of 
the projection; colours and symbols are as in Figure [T] The centre of each symbol marks the coordinates of the 
corresponding source, while the diameter of each symbol is proportional to logj^Q(DM. sin for that source. The 
greyscale is linear, ranging from to 20 rayleighs. 
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Figure 3: Best fit parameters for a planar slab of ionised gas with an exponential scale height, as a function of 
the minimum cut-off used for The top, centre and bottom panels show the scale height, i7„, the extrapolated 
mid-plane electron density, no, and the redu ced of the fit, y^, resp ectively. In the top two panels, the dot-dashed 
line shows the parameters given in Table 3 of lCordes & Lazid J2002I) . while the error bars show the 68% confidence 
intervals on each parameter In the bottom panel, the number above each datum indicates the number of points from 
Figure[T|used in the fit. The fit did not converge for a cut-off |fe|min — 10°. 
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Figure 4: The Milky Way as viewed from the North Galactic Pole, with the Galactic Centre located at the origin 
of the coordinate system and the Sun at (0 kpc, 8.5 kpc). The gre yscale represents the distribution of mid-plane 
electron density from the NE2001 model of ICordes & Lazij ("2002) (for clarity, individual high-density clumps are 
not shown). The data show the projected locations onto the Galactic {X, Y) plane of pulsars with known, reliable, 
distances. Colours and symbols are as in Figure [T] 
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Figure 5: Dispersion measure per unit distance (equivalent to the electron density averaged along the line of sight) 
as a function of height above the Galactic plane for pulsars with known, reliable, distances. Colours, symbols and 
curves are as in Figure [U 
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Table 1: Estimates of the scale height and mid-plane 
value of the volume-averaged electron density in the 
thick-disk component of the WIM. The entries have 
been sorted by increasing scale height. 



Hn (pc) 


no (cm 


Reference 


Comment 


430 - 670 


0.028 - 0.035 


Readhead & Duffett-Smith (1975) 




670+"° 
830 


0.033 ±0.002 
0.025 ±0.005 


Nordaren et al. (1992) 
Peterson & Webber r2002) 




880 ± 60 


0.019 


Taylor & Cordes (1993) 


Thick-disk component of two-component model 


oqn+250 


0.036 


Savaae et al. (1990) 




900 


0.025 


Reynolds (1991a) 


Thick-disk component of two-component model 


930 ± 130 


0.023 ±0.004 


Berkhuiisen & Muller (2008) 




950 


0.035 


Cordes & Lazio (2002) 


Thick-disk component of two-component model 


> 1000 


0.025 


Lvne et al. (1985) 


Thick-disk component of two-component model 


> 1000 


0.021 ±0.001 


BMM06 




1070- 1100 


0.018-0.020 


Gomez et al. (2001) 


Thick-disk component of two-component model 


1830ti^EJ 


0.014 ±0.001 


This paper 


Only uses data at \b\ > 40° 



Table 2: Estimates of the mid-plane volume filling fac- 
tor of the WIM (sorted by decreasing fo). 



fo 


Reference 


0.2 


Cordes et al. (1991) 


0.11 


Peterson & Webber (2002) 


0.1-0.5 


Pvnzar' (1993) 


0.1-0.2 


Reynolds (1991b. 1997) 


0.1 


Kulkarni & Heiles (1987. 1988) 


0.083 


Cox (2005) 


0.05-0.14 


Ferriere (2001) 


0.08 ±0.02 


Berkhuiisen & Muller (2008) 


0.05 ±0.01 


BMM06 


0.04 ±0.01 


This paper 



Table 3: Estimates of the mid-plane mean square elec- 
tron density of the WIM (sorted by decreasing {n^)o)- 



{n')o (cm-") 


Reference 


0.009 


Kulkai-ni& Heiles (1987) 


0.008 ±0.001 


BMM06 


0.007 


Reynolds (1990a) 


0.005 - 0.012 


Ferriere (2001) 


0.005 ±0.001 


Berkhuiisen & Muller (2008) 


0.005 ±0.001 


This paper 



